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Dark Matter via Many Copies of the Standard Model 
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We propose a cosmological scenario based on the assumption that the Standard Model possesses 
a large number of copies. It is demonstrated that baryons in the hidden copies of the standard 
model can naturally account for the dark matter. The right abundance of the hidden-sector baryons 
and the correct spectrum of density perturbations are simultaneously generated during modulated 
reheating. We show that for the natural values of inflaton coupling constants, dictated by unitarity, 
the dark-matter abundance is predicted to be proportional to the ratio of observed cosmological 
parameters: the square of the amplitude of cosmological perturbations and the baryon-to-photon 
number ratio. 



I. INTRODUCTION 

One of the major unsolved problems in cosmology is 
the understanding of the nature of dark matter ("DM") . 
The influence of its gravity can be seen not only in the lo- 
cal universe: in rotation curves of galaxies, gravitational 
lensing by them and the temperature of gas in clusters; 
but also in the form of the anisotropies in the cosmic 
microwave background. For an overview of known DM 
properties and manifestations see e.g. textbooks [3, 0|. 
All cosmological observations to date are consistent with 
the fact that the universe's energy budget comprises ap- 
proximately 4% baryonic matter, 21% dark matter and 
75% dark energy. Astrophysical data imply that if dark 
matter interacts in any other way than through gravity, 
it does so extremely weakly, leaving no observational sig- 
nature of this fact large enough to be detectable to date 
(the debate as to whether the recent PAMELA 0| and 
ATIC p| signals are a result of dark-matter annihilation 
notwithstanding) . 

Indeed, some of the most stringent constraints on the 
collider physics of supersymmetry or extra dimensions 
arise as a result of the requirement to predict the appro- 
priate dark-matter abundance. These extensions of the 
standard model ("SM") are invoked, among others, to ac- 
count for the hierarchy problem and contain within them 
natural candidates for dark matter. 

Recently, a different solution to the hierarchy problem 
was proposed in [o,0|- Namely, it was realized that in any 
effective field theory with a large number N of elementary 
particle species, the fundamental scale at which gravity 
becomes strong is lower than the Planck mass Mpj and 
is given by 
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This scale serves as the cut-off for the low-energy ef- 
fective theory. The generalization of the bound HJ to 
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quasi-de-Sitter inflationary vacua was given in Q. No- 
tice, that the scale M* also puts the upper bound on the 
mass of the species of number N. Therefore, if there 
were to exist N = 10 32 different species in our uni- 
verse, this would place the scale beyond which gravity 
is strong at 1 TeV. In this way, the strong-gravity effects 
can potentially provide a mechanism which regulates the 
ultraviolet-sensitive radiative corrections to the mass of 
the Higgs. The existence of this strong-coupling scale 
was initially proved by considering the implications of 
models with a large number of species on the consistency 
of black-hole ("BH") evaporation 0, @|. Moreover, the 
same result has been justified by considering BH entan- 
glement entropy IS and quantum information theory [§] 
and holography [jj [l3]. The possibility of the perturba- 
tive generation of the hierarchy (JJ) between the Planck 
mass and fundamental gravity scale as a result of the 
renormalization of the gravitational coupling by the loops 
of A-species was suggested even earlier [TJ LL2]. How- 
ever, in contrast to the non-perturbative BH arguments, 
the latter's one-loop mechanism had to be based on nat- 
uralness assumptions, owing to the usual limitations of 
perturbation theory. 

In this paper, we work within the above framework 
and suggest a mechanism which naturally produces the 
appropriate abundance of dark matter. We propose that 
there exist an additional, discrete permutation symmetry, 
P(N), under which the standard model is charged, and 
therefore, there exist N copies of the SM augmented by 
N copies of the inflaton. 

Our idea is that the matter we usually consider to 
be baryonic all belongs to the same species, while the 
dark matter consists of baryons of all the other copies 
of the standard model. Current astrophysical observa- 
tions constrain DM to be effectively a collisionless fluid 
that behaves as if it only interacted gravitationally. In 
our scenario, this property results from: (a) Low partial 
densities of each species except ours, which are achieved 
through the proposed inflationary and reheating mech- 
anisms; (b) The suppression of any cross-couplings be- 
tween different species compared to couplings within one 
sector, which — we demonstrate — follows directly from 
unitarity. 

The idea of dark matter's being stored in the form of 
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the several hidden copies of baryons was put forward in 
[l3| in the context of extra-dimensional theories. Since 
this mechanism was not tied to the solution of the hierar- 
chy problem, the number of copies was relatively small, 
but otherwise unconstrained. Nor was there any spe- 
cific connection predicted between the dark-matter den- 
sity and the mechanism for the generation of cosmo- 
logical perturbations. In the context of a permutation- 
symmetric standard model, this idea is also considered 
in a complementary work [1 41 ] - 

In addition, the possibility of mirror dark matter with a 
single SM mirror copy was also studied in [Tj|. However, 
as opposed to the many-copy scenario, both sectors must 
contain additional DM particles and such a possibility is 
beyond our interest. 

In other works, fltj UtI [l8| inflation with a large num- 
ber of species is also considered, however, with a different 
focus to this paper. 

In our scenario, the large number of species leads to a 
number of qualitatively novel features: 

Consistent inflation, without strong-gravity effects, 
can occur only below the cutoff M*, given by Eq. |T|). 
In this case, it is natural . I!)| to implement a modulated- 
reheating scenario (20l . l2l| to generate the appropriate 
amplitude for the power spectrum of cosmological per- 
turbations. In order to populate one of the species more 
than any other, the final stage of inflation should occur in 
one species only, so that reheating happens mostly within 
this preferred sector. However, a cross-coupling of the in- 
flaton with the other species, which naturally appears at 
least owing to gravity, ensures that some of the inflaton 
decays occur into all the other copies. The cosmological 
perturbations generated in this scenario are guaranteed 
to be adiabatic if the modulator is common to all species. 

It is important that the partial decay rates of the in- 
flaton field into different species are strongly restricted 
by unitarity considerations. As a result, dark-matter 
abundance is predicted in terms of known cosmological 
parameters — the dimensionless amplitude of the power 
spectrum of cosmological perturbations and the 

baryon-to-photon number ratio rjt,: 

fib ~ % ' (2) 

with the only additional requirement being that baryon- 
antibaryon annihilation be suppressed in the dark sec- 
tors. 

In our species, the plasma is dense and thermal equi- 
librium is achieved: it is key that the majority of baryons 
annihilate with antibaryons, creating the observed small 
baryon-to-photon number ratio. In the other species, 
given a large enough N, the individual partial densities 
are extremely low. We show that all annihilation pro- 
cesses are frozen out at all times, which ensures that all 
the dark baryons survive. 

The combination of: couplings bounded by unitar- 
ity, constraints on the modulator vacuum expectation 



value ("VEV") resulting from the observed cosmological- 
perturbation amplitude, and baryon asymmetry as given 
by the observed baryon-to-photon number ratio in our 
model predicts the correct dark-matter abundance when 
inflation occurs around the natural cut-off scale. 

The paper discusses the scenario in detail: require- 
ments for the inflationary mechanism and the mechanism 
for generating perturbations are considered in section HI] 
while the reheating process, generation of the dark mat- 
ter and the subsequent evolution of the baryonic and dark 
sectors are laid out in section Hill 



II. INFLATION AND PERTURBATIONS 

A. Inflationary Background 

We consider a theory with a large number, N, of stan- 
dard model copies, which we label by an index j — 
1,2,...N. For maximal predictivity and simplicity, we 
assume that each of the copy sectors is an exact replica of 
the SM, and all such sectors are related by an exact per- 
mutation symmetry P(N). Motivated by the suggested 
solution of the hierarchy problem 0,0], we keep in mind 
that the number of species can be as large as N ~ 10 . 
However, we show that the scenario proposed here also 
works for a range of values of N. Therefore, for gener- 
ality, we will keep TV as a free parameter in most of the 
below. 

It is widely accepted (see textbook discussion e.g. 
P, 0, that the history of the early universe starts 
with a period of inflation driven by an effective scalar 
degree of freedom $, which has a sufficiently flat poten- 
tial that allows the field to slow roll. Taking in mind the 
benefits of inflation we shall not question this paradigm, 
and supplement the SM fields by the inflaton 1 . Our N- 
species framework allows us to endow each j-th SM copy 
with its own inflaton field $ J ; enumerated by the same 
index j. 

At the level of the lagrangian, there is a full permuta- 
tion symmetry between the SM copies and the associated 
inflatons. It is clear from our experience of the world that 
this symmetry must be broken by the cosmological back- 
ground since we are coupled to a fraction larger than one 
in N of the matter in the universe. We posit that the seed 
of this symmetry breaking was prepared by inflation. All 
we require is that the evolution of the last few e-foldings 
of inflation be driven by a single inflaton species. Such 
a situation is much more natural than one may naively 
think: 

The reason is the bound HJ applied to the inflationary 
Hubble parameter, which tells us that, unless the indi- 
vidual inflatons give small (<C M») contributions to the 



It would be interesting to consider whether this role could be 
played by the Higgs field. 
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Hubble parameter, inflation can be driven at most by 
a few inflatons simultaneously. To illustrate this point, 
consider a simple case — the potential, 



(3) 



where equality of all the individual potentials arises ow- 
ing to the permutation symmetry. We assume that 
each individual sector is capable of driving inflation in 
the absence of contributions from the others, just as in 
minimal inflationary scenarios: for the right choice of 
inflaton VEV, the would-be Hubble parameters Hf = 
V($j)/3Mpj are not much below M* and the potentials 
have appropriate slow- roll parameters. 

In such a case, there would exist an inflationary 
state where inflation is being driven by all N infla- 
tons, giving an effective Hubble parameter of — 
Ej^( $ j)/ 3M pi ~ NM *- This is far in excess of the 
bound ([I]) and therefore inconsistent. The system should 
contain a built-in mechanism to prevent such a collective 
inflationary state: we can add cross-couplings between 
the species to achieve this, for example, 



v = E v ^ + E K*. 



+ Jfi E + hi S her terms - ( 4 ) 



In the case of a cross-species coupling, there are new con- 
straints that arise as a result of the requirement for the 
perturbative unitarity of the theory. For example, con- 
sider the simplest cross-coupling A 2 2 < I ) f Having such a 
coupling allows us to consider annihilation diagrams for a 
particular species i to a species j mediated by loops of all 
other species (see Fig. [I]). Since the species index on each 
loop is free, such diagrams will involve summations over 
all the species, and therefore will be relatively enhanced 
with respect to one-specie theories. Diagrams with n con- 
secutive loops are possible and will scale as A22(A22-V) n . 
This indicates that we must constrain A22 < N^ 1 in or- 
der to prevent the amplitudes from diverging. 

Similar bounds follow from many other considerations 
including the renormalization of the kinetic terms of the 
species due to loops of other species. This scaling is 
also fully compatible with the typical values of the cross- 
couplings induced via the exchange of virtual black holes 
: the dominant contribution comes from black holes of 
mass M* and each exchange by a virtual black hole be- 
tween the different species requires a suppression by 1 /N 
in the amplitude. 

Let us assume that all the couplings are positive- 
definite. In that case, despite the required suppression 
of the cross couplings, putting one of the inflatons on the 
inflationary slope at some VEV, say > Mpj, forces all 
the other inflatons to remain at zero VEV or, if they had 
a non-zero VEV before, to occupy zero VEV during the 
time of order H . Indeed, the effective mass term for 




n-loops 

Figure 1: Introducing a cross-species coupling to the infla- 
ton potential involves diagrams which couple species through 
loops. Since the species index on the loops is free, these am- 
plitudes are relatively enhanced by a factor N for every loop, 
M ~ A22(A22A r ) n . Perturbative unitarity of the theory there- 
fore requires that the coupling A22 be suppressed by at least 
a factor N . 



all inflatons with index j ^ k in the background of is 
given by 
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This effective mass is at least of order M* , and thus all 
the other inflaton fields strongly violate the slow-roll con- 
dition: the effective friction cannot prevent them from 
fast roll towards the minimum. At the same time, quan- 
tum perturbations for the degrees of freedom with mass 
M<s> jllk > H* are not enhanced, leaving no way to pro- 
duce a VEV for them. 

Moreover, non-perturbative black-hole arguments can 
be used to derive a bound in (quasi) de-Sitter spaces on 
the masses of particles which have N species [3], 



m < 



M P1 



N 



(6) 



Any inflationary vacuum in which this bound is not sat- 
isfied must destabilize and cannot support slow-roll. If 
the coupling A22 exceeded 1/N, the inflationary vacuum 
with $fe ^> Mpi would result in the presence of order 
N species with masses 3> Mpi/s/N, violating the bound 
above. This argument is fully non-perturbative and is 
independent of the unitarity considerations, but results 
in the same constraint on the cross-species couplings as 
those we found through unitarity argument above. 

The system of cross-couplings also eliminates the unde- 
sired possibility of inflating collectively with all inflatons 
at small VEVs, but with £\ $| — Mp,. Such a situ- 
ation also provides sufficient inflation for at least some 
forms of the potential Eq. |J5J, but reheats to all the 
species equally. With the cross-couplings present, the 
cross-terms add up in such a way that the mass of all the 
inflatons is high and no inflation can occur. 

Overall, despite the symmetry of the potential, the 
presence of the cross couplings ensures that inflation can 
only proceed in one of the inflaton species but not several 
simultaneously. 
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We can further generalize the cross-coupling terms in 
Eq. 01): 



Here, the first term is the sum over the individual inflaton 
potentials, while the second term is the sum over the 
cross couplings, meaning that in each term at least two 
or more distinct species participate (that is n > 1). In the 
presence of individual quantum numbers of the inflatons, 
the powers ctk are constrained by symmetries. 

As before, unitarity arguments constrain these multi- 
specie cross-couplings to scale with N as X ai ,...a n J$ 
iV 1- ™, i.e. couplings between n different species need to 
be suppressed to at least a factor N 1 ~ n since, for each 
species, a loop can be formed and summed over. 

The particular choice of potentials discussed here and 
the chosen cross-couplings serve only an illustrative pur- 
pose, capturing the essence of the mechanism through 
which the fundamental theory can prevent collective in- 
flation in order to obey the bound on the Hubble param- 
eter. From here on, we shall keep the form of the individ- 
ual inflaton potentials fully generic. For the purpose of 
model building, notice that there is no inconsistency in 
having inflaton VEVs significantly above M», as long as 
the value of the Hubble parameter does not exceed the 
fundamental scale M* 2 . 



B. Perturbations 

The main benefit of inflation is a natural quantum 
mechanical generation [23] of primordial cosmological 
perturbations which are necessary to seed the inhomo- 
geneities observed in the universe. For a modern peda- 
gogical treatment of this mechanism and historical details 
see e.g. 0,0, El- 

The generated inhomogeneities are usually character- 
ized by the dimensionless power spectrum A^(fc) of cur- 
vature perturbations 1Z. In models of inflation involving 
a single scalar field with a canonical kinetic term, this 
power-spectrum is given by: 
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parameter, and we have used the reduced Planck mass 
Afpi = (8nG N )~ 1/2 . The function A^(fc) has usually 
only a slow logarithmic dependence on the scale k which 
can be approximated in the narrow observable range by 
a power law A^(fc) oc k ns ~ 1 . One of the spectacular suc- 
cesses of observational cosmology is the measurement by 
the Wilkinson Microwave Anisotropy Probe ("WMAP") 
[25[ of the spectral index, ns = 0.960 ± 0.013, and the 
value of the power spectrum, A^ = (2.445 ± 0.096) x 
10 -9 , at the scale k — 0.002 Mpc -1 . This measurement 
is in an excellent agreement with predictions of the sim- 
plest inflationary models. 

For a larger number of inflaton species, the result of 
Eq. (HI) is not modified [IB, H3|, see also [2, page 506]. 
Thus, in chaotic inflation, in the case of many species, 
where the Hubble constant is bounded by the fundamen- 
tal scale |H), one obtains 
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where is the Hubble parameter at the time of horizon 
exit of a particular mode k, e = —H/H 2 is a slow-roll 



From this inequality, it follows that only models with 
'moderate' values oi N (N < 10 10 ) can avoid extremely 
small values of the slow-roll parameter e and potential 
tensions with the measurement of ns quoted above. This 
bound for the number of species in inflation can be ob- 
tained in different ways see e.g. [7|, ll6j, [28], l29j. In theo- 
ries with these moderate numbers of species, cosmological 
perturbations can be generated in the standard manner 
and no alternative mechanisms are required. However, 
for larger values of N, inflation is forced to proceed at 
a lower scale and the amplitude of density perturbations 
generated in standard chaotic inflation becomes unac- 
ceptably small. 

The issue of generating density perturbations in the 
presence of a large number of species is studied in [l9| . 
A number of successful scenarios exist, but the mecha- 
nism most suitable for our purposes is modulated reheat- 
ing (20L Hl|, which enables the generation of sufficient, 
red-tilted cosmological perturbations for the small val- 
ues of -£/* required by the bound Eq. UJ . In this mech- 
anism, one introduces a new scalar field x, the modu- 
lator, which controls the decay rate, T, of the inflaton 
field into SM particles. This modulator \ is assumed to 
be light, M x <C -£/*, so that substantial quantum fluctua- 
tions 8\ ~ are generated during inflation. These fluc- 
tuations become frozen (classical) as they leave the hori- 
zon. Moreover, it is assumed that x does not contribute 
significant energy density during inflation and does not 
affect inflationary dynamics: 



2 The question of whether the inflaton VEV can exceed Mpi is 
separate and more profound. The impossibility of this can be 
proven in certain cases and results from the conflict with black- 
hole thought experiments that one can perform in de-Sitter space 
Q. In the context of this work, we will not worry about this 
issue, since inflation models with sub-Planckian field VEVs can 
be constructed. 
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As usual, the inflationary epoch is followed by the oscil- 
latory dust-like stage. During this stage the Hubble con- 
stant decreases quickly, e > 1, and eventually becomes 
less than the decay rate T. Around that time reheat- 
ing occurs: inflaton decay becomes efficient enough to 



populate the universe with SM particles. Without loss 
of generality, we assume that during reheating the infla- 
ton(s) oscillate about $j = 0. The local VEV (x) con- 
trols the decay rate of the inflatons, T, and therefore con- 
trols the position of the reheating surface. Perturbations 
in the field value of the modulator lead to perturbations 
in the reheating surface. Assuming that the universe was 
dust-dominated prior to reheating, radiation-dominated 
afterwards and that reheating occurs instantaneously at 
H ~ r, one obtains [13, HH for the curvature perturba- 
tion in the radiation fluid 
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where the final equality results from the fact that the 
couplings which provide the dominant decay channels are 
linear in the modulator field. We have also assumed that 
the amplitude of perturbations in the inflatons is negli- 
gible compared to that imparted to the modulator. 

We require that the modulator \ be common for all 
species. Thus all the perturbations in the decay products 
of the inflaton are generated by a single source: pertur- 
bations in the modulator field. This ensures that, even 
if the decay products are not thermally coupled after re- 
heating, their curvature perturbations will be identical 
and will not contain an isocurvature mode. The relative 
contributions to the energy density in different particle 
types will be the same everywhere, only the time of re- 
heating will have changed. 

The modulation of the reheating process can be guar- 
anteed by invoking a symmetry which ensures that all 
inflatons couple to the light SM fields through the mod- 
ulator x- The simplest possibility is to require a Z2 sym- 
metry that we shall refer to as inflaton parity. We shall 
assume that this symmetry is common for all the copies. 
Under this symmetry, all the inflatons, as well as the 
modulator, change sign, 
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while all the usual SM fields and their replicas are left 
invariant. This fact ensures that the decay rate of all the 
inflatons is modulated by the fluctuations of the expec- 
tation value of x- 

The dimensionless curvature power spectrum resulting 
from this reheating mechanism is given by 
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For the correct amplitude of the curvature perturbations, 
the vacuum expectation value for the modulator during 
inflation is then required to be 



III. REHEATING AND FREEZE OUT 

A. Inflaton Decay 

The lowest-order invariant couplings between the in- 
flatons and the SM species (call them Qj) can be param- 
eterized in the following way 



gQjQj) , 



(14) 



where the two coupling constants g and g are respon- 
sible for the diagonal and non-diagonal decays respec- 
tively. Notice that, owing to the underlying structure 
of the theory, these effective coupling constants may de- 
pend on some VEVs. For example, if Q are quarks and 
the inflaton is a SM gauge singlet, the coupling constants 
may depend either on the Higgs VEV or on powers of 
the four-momentum. The convenience of the above pa- 
rameterization is that such details are unimportant: we 
express everything in terms of the cutoff and the effective 
dimensionless couplings which must satisfy the unitarity 
constraints, irrespective of their substructure. 

In our analysis, we shall assume that the values of these 
coupling constants are the most natural (the largest pos- 
sible) compatible with the requirement of unitarity at the 
energies below the cutoff scale. We discuss these unitar- 
ity bounds below. 

As was already explained in the previous section, the 
general iV-scaling of the inter-species couplings is such 
that an operator involving n different types of species 
must be suppressed by l/N 71 ^ 1 . The most stringent uni- 
tarity bound on g comes from this consideration. Even 
though this operator is diagonal in the species index, it 
can generate inter-species interactions via the exchange 
of a virtual x- Such a process is shown in Fig. EJa). 
At characteristic energies E < M x , we can integrate the 
modulator out to obtain an effective six-point interac- 
tion between the two types of species, with the vertex 
g 2 /M 2 M 2 . Assuming that the mass parameters are of 
roughly the same order, we get the following constraint, 
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The bound on the non-diagonal coupling, g, is at least 
1/N, since it involves species with two different labels. 
However, a more stringent constraint can be obtained 
since the modulated-reheating mechanism requires the 
presence of a background field VEV, (x) • On such a back- 
ground, we can reduce the g vertex of Eq. 1|14[1 to a three- 
point interaction. The unitarity bound then arises from 
the diagram involving the exchange of virtual-inflaton 
species, shown in Fig. (Hb), with the effective coupling 
Ng 2 (x) /M%M%, where the factor TV results from the 
summation over all the species of virtual inflaton. The 
constraint on the coupling then becomes 
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set by the ratio of the inflaton decay rates into the two 
sectors 
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(P) 



Figure 2: Processes which lead to strongest constraints on 
the coupling constants g and g arising as a result of interac- 
tion terms Eq. Ill4ll . Bounds are governed by the requirement 
for perturbative unitarity. (a) Exchange of a virtual %■ 00 
Exchange of a virtual inflaton $fc; total amplitude involves a 
sum over all inflatons k. 



Using (fl~3|) . this constraint can be recast as 
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One can also constrain the maximal inflaton VEV dur- 
ing inflation. Namely, unitarity arguments are most re- 
strictive for the rate of the process Qi + Qi Qj + Qj 
to any j, occurring in the presence of For the char- 

acteristic energies E, the rate is bounded in such a way 
as to satisfy T(E) < E. Demanding that unitarity hold 
up to the cutoff scale and using the constraints on g and 
g, bound l|15p and bound (fT7|) respectively, one obtains 




100. 



(18) 



where we have used Eq. (jT3J) . Or, equivalently, the total 
energy density in the dark sector right after the reheating 
is 



pj 



A 2 
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(20) 



and is strongly suppressed compared to the energy den- 
sity in our sector. Note that this result does not depend 
on the number of species N . Since interactions between 
the sectors are suppressed, the other species are never 
in thermal equilibrium with our sector. This means that 
they are not included in the naive tally of energy density 
in radiation and therefore contribute an energy density 
that is unaccounted for during Big-Bang nucleosynthesis. 
The smallness of Eq. lj20|) guarantees that the nucleosyn- 
thesis bound on the amount of extra radiation in the 
Universe [30l is satisfied. 



B. Reheating in Our Sector 

Let us now estimate whether we can reheat the uni- 
verse to a sufficiently high temperatures Trh to synthe- 
size the light elements in successful Big-Bang Nucleosyn- 
thesis ("BBN"). The reheating takes place once the decay 
rate of the inflaton T^> 1 ~ H. The branching ratios to 
species other than j — 1 are subleading, so we can use 
the approximation that the total decay rate is equal to 
the decay rate to our species, 
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is the decay rate of the inflaton: $i — ► Q± + Q\, There- 
fore, right before reheating, the energy density stored in 
the inflaton $i is 
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Note that this bound is not restrictive at all. 

As described in section |Ul the decay products after re- 
heating all come from a single inflaton species, $i. The 
coupling g is relatively large, Eq. lfT5]). and controls the 
diagonal decay into matter species Q\. However, parti- 
cles in the other copies are also produced by virtue of the 
much smaller off-diagonal coupling g, Eq. (fT7|) . The in- 
dividual energy densities in all the other sectors are neg- 
ligible relative to our SM sector; however, the integrated 
energy density is significant. The ratio of the energy den- 
sities between the SM and any other given j-th copy is 



where we have used Eq. (|T5j) and Eq. (|13|) . Then, within 
time ~ rr , the energy of inflaton oscillations gets con- 
verted into the energy of the species. During the decay, 
Qi are created with the energy of order E ~ M$ x . This 
is true for all the species. The number density in our 
sector, n\, is large and allows for rapid thermalization. 

Let us estimate the thermalization rate. This is the 
same as the interaction rate of SM particles: 
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Since, from (|22|) the Hubble parameter at reheating is, 



h rh ~ ri ~ m* i -f 
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(24) 



we obtain the result that the number of scatterings per 
particle Qi during one Hubble time is 
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Using ((TJ this expression above can be rewritten in the 
form 



r T 



-Hr.h A tc \M $1 
In our inflationary scenario, H* > , hence 



r T a 2 

> TT" > 1) 



#RH ~ A| 



(26) 



(27) 



The reason for the final domination of the dark-matter 
sector over the baryons lies in the difference in their be- 
havior resulting from the huge hierarchy in the partial 
energy densities of the individual sectors. 

The thermal history for our SM sector is unaltered: 
The high density of SM baryons and anti-baryons en- 
sures that they remain thermally coupled, converting 
the baryons and antibaryons into photons extremely effi- 
ciently. Only as a result of some primordial baryon/anti- 
baryon asymmetry does a significant number of baryons 
survive. Nonetheless, these represent a tiny fraction of 
the original number, which is known to be 



Vb 



^=2.7-10- 



n h h 2 =6.2- 10" 



(31) 



This number is also known to not have been significantly 
different as early as at the time of nucleosynthesis 0|. 
Therefore, for every baryon in our species that exists to- 
day, approximately r/^ 1 baryons in our copy of the SM 
were in existence after reheating. 



and the thermalization happens well within a Hubble 
time. Therefore, the reheating temperature can be es- 
timated using 



RH 



t rh r 9*<X) M 9l 



M P1 " 1 Ml 
which combined with Eq. lfl3|) yields 
9 (H, 



(28) 



(29) 



Further, we apply the constraint Eq. I|15p along with the 
bounds on the Hubble parameter (1), < M*, and on 
the inflaton mass Eq. |(6|), < M*, and obtain the 

following upper bound for the reheating temperature 



RH 



< Mli N -3/4 



10 5 M P ,iV~ 3/4 . 



(30) 



For N ~ 10 32 the reheating temperature is around or less 
than 1 GeV and it is only higher for smaller N. In partic- 
ular, for N ~ 10 20 we have T RH ~ A/* - M $1 - 10 9 GeV. 
It is worth mentioning that Trh ~ 10 5 M t ,N~ 1 / 4: , there- 
fore for all N > 10 20 the reheating temperature is much 
less than the assumed mass of the inflaton < M« . 
The number of species N clearly cannot exceed 10 32 by 
more than a few orders of magnitude, since otherwise 
the strong-gravity effects would have already been seen 
by accelerators. Hence, in our scenario Trh > 1 GeV 
for all interesting N, which ensures that BBN proceeds 
unaltered. 

Thus, shortly after reheating, the universe is left with 
our sector at temperature Trh, while the other sectors — 
with much lower particle densities — remain unthermal- 
ized and consist of a collection of particles with energies 
- M 4l . 



C. Dark Sector Freeze-Out 

For the baryons not in our species, the story is very 
different: given a large enough N, their partial densities 
following reheating are extremely low. In our model, we 
require that they be so rare that all annihilation processes 
are always frozen out and the other species never reach 
thermal equilibrium. We define freeze-out as the time at 
which the annihilation rate for the j-th copy, T*, falls 
below the Hubble parameter. 

Let us consider the annihilation of hidden protons and 
antiprotons in the species j at the time when the the tem- 
perature of the particles in our sector is T. The number 
of annihilations per particle during one Hubble time is 
given by 



H 



UjUV 

H 



(32) 



where rij is the baryon number density of the dark species 
in sector j, a is the annihilation cross-section and v is the 
typical speed of the particle. The density rij is conserved 
in a comoving volume, provided that annihilations be 
suppressed. We can, therefore, express nj in terms of the 
density of our baryons as 



rij(T) = m 



Pi \Trh 



(33) 



where m, Pj and p\ are evaluated before thermalization 
occurs. In particular, the number density in our species 
can be estimated using Eq. (|22| and Eq. ([28)1 



P>i> 



M $1 \ M *i 



RH 



rp3 

RH' 



(34) 
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Thus, using Eq. (|T9|) . we obtain the j baryon number 
density as the function of the temperature in our ther- 
malized sector 



nj(T) = 



RH 



A 2 
N 



(35) 



There are three primary channels through which mass 
stored in protons can be reduced: (a) annihilation 
through a virtual photon to charged particles (electrons, 
positrons, pions, etc): jj — ► e~e + , (b) annihilation to 
two photons: jj — ► 27, (c) formation of a bound state 
through a radiative-capture process followed by annihila- 
tion: jj — > bound. There are also three regimes which we 
will need to consider for the cross-sections: relativistic, 
non-relativistic and ultra non-relativistic or Sommerfeld 
enhanced [3l]]. As we have already mentioned, at pre- 
heating, the characteristic energy of a proton in the j-th 
copy of SM is £(Trh) ~ M$ 1 ^ m p, where m p is the 
mass of the proton. Thus, copies of the protons are ini- 
tially highly relativistic. 

During the expansion of the universe, the momentum 
of the protons in the j-th copy of the SM redshifts, p oc 
a -1 . This implies that, while they are relativistic, the 
characteristic energy redshifts as E(T) ~ M^T/Trh, 
Eventually the kinetic energy becomes comparable to the 
rest mass of the proton and it becomes non-relativistic; 
this occurs when the temperature in our sector is 



Ah 



T, 



RH 




(36) 



At temperatures below T cr , the momentum continues to 
redshift which leads to a redshifting of velocity, v ~ 
T/T CI . Finally, when the speed of the dark-sector protons 
falls beneath the fine-structure constant, v < a, Som- 
merfeld enhancement becomes important for the cross- 
sections [1H. The temperature in our sector at this tran- 
sition is 



21 



Som 



(37) 



We will now consider the freeze-out condition for each 
of the three annihilation channels in turn. Using H ~ 
T 2 /M Pl and Eq. $35$) in Eq. (32]) yields 



H 



^) "PIT. 



(38) 



We have presented the parametric forms of the cross- 
sections av in Table HI Using these in Eq. |38|) allows 
us to investigate the evolution of the freeze-out condition 
for each of the channels. 

a. Annihilation to charged particles jj — > ee + . In 
the relativistic regime, T > T cr , we can estimate 
a ~ a 2 /E 2 , by redshifting the energy we obtain av ~ 
Substituting this into Eq. lj38j) yields 



a 2 T 2 J {M^Tf 



(jj->ee+) 

H 



2 A K ( Tl 



RH 



pi 



T 



(39) 



Regime 



Relativistic 



Non- Relativistic 



Sommerfeld 



Charged 



Photons 



Capture 



a 
E 2 



m p E 



In 



2E 



a 
m 2 



( 9L\ a 2 ( a\ 
\ v J m 2 \ v ) 



a 2 /a\ 4 
m 2 \ v J 

4 (-) 



Table I: Parametric forms of the cross-sections, av, for three 
potential annihilation channels for protons and antiprotons: 
(a) annihilation to light charged particles through virtual pho- 
ton [32| ; (b) annihilation to two photons [32]; and (c) forma- 
tion of a bound state through radiative capture [331 ] . We 
present three limits of the expressions: highly relativistic, 
valid for temperatures in our sector T cr < T < Trh, non- 
relativistic, valid for Ts om < T < T cr and Sommerfeld en- 
hanced [HI], valid for T < T Som . 



Here it is important to note that as a result of the 
redshifting of the relativistic E(T) the above function 
grows as the universe cools down. When the temper- 
ature in our sector reaches T cr the dark baryons en- 
ter the non-relativistic regime. There, according to Ta- 
ble [U av = const and consequently Ffjj_ te6 +)/ H de- 
creases as T. Finally, after the universe cools to Tsomi the 
Sommerfeld-enhanced regime starts, and the amount of 
annihilations per particle per cosmic time, T^j^^+^/H 
approaches a constant owing to an additional power of v 
in the denominator of the cross-section. Thus, the max- 
imum of Trjj^ ee +^/H is reached when the dark baryons 
are on the verge of becoming non-relativistic, at T = T„. 
If this annihilation channel for the dark baryons is frozen 
out at this maximum, then this annihilation reaction is 
frozen out at all times. 

It is convenient to rewrite the expression lj39"ll using the 
formula JH and Eq. (fT5)l 



T (n- 



-) a 2 ( M] 



pi 



H 



N 2 



RH 



T 



(40) 



Substituting the critical temperature in the freeze-out 
condition yields a bound on the lowest number of species 
necessary to ensure that annihilations to light leptons in 
the dark sector be frozen out, 



N > 



,Mpi M* 
m p M$ 



2/3 



10 a 



(41) 



b. Annihilation to two photons jj — ► 27. Analo- 
gously to the previous case, for the relativistic regime, 
we use the cross-section from Table [J and the redshift 
law for the energy, allowing us to write Eq. lj38j) as 



r 



( jj->27) 

H 



a 2 A|, Mpi 
N m n 



Xrh 



In 



2M $1 T 
TOpTRH 



Thus in the relativistic regime T^j^ 2 -y)/H decreases log- 
arithmically as the universe cools. However, the non- 
relativistic regime is identical to that of process (a), i.e. 
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T(jj^ 2 -y)/H first decreases as T and then becomes con- 
stant, once Sommerfeld enhancement is active. There- 
fore, T(jj^ 2 y)/ H is maximal at Irh and consequently N 
will be most strongly constrained at reheating where 



(jj->2 7 ) 



RH 



TV 3 / 2 \M<s,J V m p 



M P1 \ , / 2M$, M P1 
In 



M* m p ViV 



Requiring that this process be frozen out at reheating, 
constrains the number of species to be 



N > 10 1 



(42) 



c. Formation of bound state. Once the protons are 
non-relativistic, it is possible for them to be captured 
into a bound state. For the purpose of this estimate we 
have used the asymptotics of the cross-sections obtained 
in [33]. Once the bound state is formed, the decay rate of 
the proton-antiproton pair is extremely high, ~ m p a 5 ~ 
10 10 GeV [H. 

It is clear from the cross-sections presented in Ta- 
ble U that, for v > a, the cross-section for capture is 
highly suppressed compared to the annihilation channels 
(a) and (b) and therefore the formation of the bound 
states does not occur provided those channels be frozen 
out. Once the universe cools further and reaches v < a, 
r '( j] abound) /H becomes constant of the same order as 
that for the other channels. Therefore there are no new 
bounds on N arising from this channel. 

d. Annihilation in halos. Another important factor 
to consider is the annihilation of dark protons in galactic 
halos. Once the halos are formed as a result of non-linear 
collapse, the density inside them is much higher than the 
cosmological energy density. This makes it possible for 
baryons to interact much more efficiently. At T < Ts om , 
the cross-sections for all the channels discussed above 
are parametrically equal (see Tabled, therefore all these 
processes will place the same constraint on the number 
of species. 

We are going to estimate an upper bound for the rate 
of the annihilation processes. A halo can be defined as a 
region with average density p^ain > 200p m , where p m is 
the background matter density [34|. We will assume that 
from the redshift at which the halos have first formed, 
z form ; until today, all the mass of the universe is stored 
in halos. In such a circumstance, we can write down a 
redshift-dependent annihilation rate 



halo 



Pha,\o(z)a 3 
Nm 3 v(z) 



(43) 



For the purpose of this upper bound, we will also assume 
that, throughout the period under consideration, the ve- 
locity of protons does not increase as a result of virial- 
ization, but remains constant at the value it had before 
the dark protons fell into the halos, v ~ aT(z{ oxm )/Ts om . 
With these assumptions, using Eq. (|37| we obtain the 
result that 



r 



halo 



(z) 



200p (l + z) 3 a 3 
i 2 p T(z (mm )A n N^ 



(44) 



Since T^i^z) represents the rate per particle of interac- 
tions at redshift z, we can ensure that the annihilation 
in halos remain insignificant by requiring that the time 
integral of Eq. f44|) over the period of the existence of the 
halos be small, say 



t(z=0) 



t(z=Zform) 



halo 



{z(t))dt < 0.01. 



(45) 



Using the conservative value of Zf or m ~ 10, the Hubble 
parameter today Hq ~ 10~ 34 eV and To ~ 10 -4 eV we 
obtain the bound 



N > 10 



li 



(46) 



Putting all of the above results together, we find that 
the strongest bound on N comes from the annihilation to 
two photons at the time of reheating. Hence our scenario 
works for the number of species in the window: 



10 11 < N < 10 



:!2 



(47) 



Nonetheless, provided that some additional mechanism 
which suppresses the interactions between the dark 
baryons existed, the annihilations could be frozen out for 
much smaller N, still giving the appropriate dark-matter 
abundance. 

Here, a cautionary note is necessary. In the anal- 
ysis above we have naively assumed that the inflaton 
decays practically directly to a small number of dark 
(anti)protons. This is not unreasonable for Mq, 1 ~ M* ~ 
1 TeV or, equivalently, N ~ 10 32 . However, for a smaller 
number of species N, the fundamental scale is far above 
TeV and new physics is required once again to solve the 
hierarchy problem. This new physics could make our as- 
sumptions invalid. Another potential deviation from our 
setup could be caused by the fact that decays of ultra- 
massive inflatons may result in the production of a very 
large number of protons per initial inflaton particle. Each 
dark proton would then carry much less energy than the 
inflaton mass, contrary to the assumptions that we have 
made. This would enhance the annihilation rates and 
therefore make the lower bound on N higher. In par- 
ticular, if the energy per proton after inflaton decay is 
M^ 1 /fj,, then the expressions for the number of annihila- 
tions per cosmic time, e.g. the right-hand side of Eq. (j39"l) . 
are enhanced by /i 2 . For example for N = 10 32 and 
M $1 ~ M*, we have T {fj / H (T RH ) ~ 10~ 33 /i 2 while 
fi cannot exceed 1000. Thus even a very high multiplicity 
would not change our conclusions and can at most affect 
the lower bound on N. 

It may seem surprising that dark matter could be bary- 
onic: in the presence of large densities of baryons, the in- 
teraction rates are large and they are able to cool and fall 
into the center of halos to form stars. However, given a 
large enough number of copies of the SM, N, each species 
becomes so dilute that the probability for any baryon to 
meet another from its own species becomes extremely 
small. Since the nucleon density today is approximately 



10 



n° = 3 ■ 10~ 7 cm' 3 , with N = 10 32 we would have 
a partial density in any other species of approximately 
~ 10~ 38 cm -3 . Unitarity bounds cross-specie cou- 
plings to be extremely small and the baryons become 
nearly collisionless, interacting effectively only through 
gravity. 

The limits on the mass for such dark-matter candidates 
are weak: constraints arising from the Lyman-a forest 
limit the mass of a sterile neutrino to be m v > 10 keV 
fUjllHl- In our proposal, the majority of the mass of the 
universe is in protons and anti-protons in other species, 
with their mass m p ~ 1 GeV easily satisfying the current 
constraints. 

Finally, we can calculate the contribution of such dark 
matter to the total energy budget of the Universe today. 
Taking the result Eq. ([20"1) for the ratio of the density 
in our sector to the dark-matter density and adjusting 
it for the fact that only a fraction given in Eq. ([3ljl of 
baryons in our sector will fail to annihilate to radiation, 
we obtain, 

This is rather remarkable: taking the inflationary Hub- 
ble parameter and inflaton mass to their most natural 
value, around the cutoff, together with the saturation of 
the unitarity constrains on coupling constants gives the 
correct dark-mater abundance. 



IV. CONCLUSION 



This field mediates couplings between the inflaton and 
both our sector and the other species. The observed am- 
plitude of cosmological density perturbations determines 
the value of the VEV of the modulator. Combining this 
with the above-mentioned unitarity constraints produces 
a natural hierarchy between the inflaton couplings to our 
sector and to the other species. The amount of matter 
produced in each species is determined solely by these 
considerations. 

In our species, the density is large and therefore 
baryon-antibaryon annihilation proceeds unimpeded. 
Provided that the number of species be sufficiently large, 
for example appropriate to account for the hierarchy 
problem, the annihilation in the other species is always 
frozen out. 

The end result is that the final ratio of dark-matter 
to baryon densities is equal to the ratio of the amplitude 
of cosmological perturbations and the baryon-to-photon 
ratio, given that inflation occurred at the natural scale 
of the cut off of the theory. 

It would be very interesting to find a mechanism which 
would be able to relate baryon asymmetry to the am- 
plitude of cosmological perturbations, or other models 
which could connect these to the DM abundance. In 
addition, the possibility exists that the dark matter in 
our model may, in fact, be millicharged [37j and couple 
weakly to the photon in our sector. This opens an avenue 
for potential direct or indirect detection. Finally, it is in- 
teresting to note that, for a sufficiently large N, electrons 
also survive and could eventually form dark atoms with 
the protons. 



In this paper we have shown that, in models with 
a large number of species, we can naturally reproduce 
the observed dark-matter abundance. The parameters 
of such models are constrained by the requirement for 
perturbative unitarity: the effective coupling between 
species must be suppressed in such a way that ampli- 
tudes with loops of species do not diverge. 

We have introduced a permutation symmetry, which 
results in our Universe's containing N copies of the stan- 
dard model. Inflation occurs through the inflaton be- 
longing to our copy, while perturbations and reheating 
are driven by a modulator field common to all species. 
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